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Abstract 

We confront the recently proposed exact S-matrices for AdSz/CFT2 with direct 
worldsheet calculations. Utilizing the BMN and Near Flat Space (NFS) expansions 
for strings on AdS 3 x S 3 x S 3 x S 1 and AdS 3 xS 3 xT 4 we compute both tree-level and 
one-loop scattering amplitudes. Up to some minor issues we find nice agreement in 
the tree-level sector. At the one-loop level however we find that certain non-zero tree- 
level processes, which are not visible in the exact solution, contribute, via the optical 
theorem, and give an apparent mismatch for certain amplitudes. Furthermore we find 
that a proposed one-loop modification of the dressing phase correctly reproduces the 
worldsheet calculation while the standard Hernandez-Lopez phase does not. We also 
compute several massless to massless processes. 
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1 Introduction 



The AdS 3 /CFr 2 -correspondence, where the gravity side consists of supcrstrings on either AdS 3 x 
S 3 x T 4 or AdS 3 xS 3 xS 3 xS 1 supported by pure Ramond-Ramond (RR) flux and preserving 16 
supercharges, has recently been seen to have an interesting integrable structure 1 - 12] 1 similar 
to that of strings in AdS 5 x S 5 and AdS^ x (DP 3 . In the more general S 3 x S 1 case very little 



is known about the dual CFT 14 15 but one can still learn a lot by analyzing the integrable 



structure of the theory. On the string side one can use the standard approach of a supercoset 
model description, to formulate sets of quantum Bethe equations [T 16-20] (see also [H] for 
a discussion of integrability without fixing the kappa-symmetry of the Green-Schwarz string). 
Alternatively one can start from the superisometry group of the background and construct a 
spin-chain from which one can read off the S-matrix (which in turns is diagonalized by the Bcthc 
ansatz equations) [2l](2^p ppl] . The relevant algebra for AdS 3 x S 3 x S 3 x S 1 is d(2, 1; a) 2 (this 
corresponds to a subgroup of the superisometry group of the background D(2, 1; a) 2 x U(l)). 
The parameter < a < 1 sets the ratio of the two S 3 radii [l]. The T 4 corresponds to a = 0, 1 
and in this case the relevant algebra becomes psu(l, 1|2) 2 = d(2, 1; 0) 2 = d(2, 1; l) 2 . 

In this paper we continue the investigation initiated in [4j[5] in order to confront the proposed 
S-matrices of [7] and [8 ,11 with explicit worldsheet calculations. We do this both at the tree 
and one-loop level utilizing the BMN and NFS expansion of the string in AdS 3 x S 3 x S 3 x S 1 

24 , see also [25]. The worldsheet excitations consist of four pairs of 
(1, a, 1 — a, 0) and at tree-level we compute all possible processes for 



23 



and AdS 3 x S 3 x T 4 
modes with masses m 

two light (m — a, 1 — a) incoming bosons. For the case of AdS 3 x S 3 x T 4 , where we have a four 
fold degeneracy of heavy (m = 1) and massless modes, we also perform a one-loop computation 
in the near flat space (NFS) limit. Upon natural choices of free parameters and some minor 
adjustments we find that the tree-level sector of the S-matrix found in [8 11 reproduces our 



1 For the case of both NSNS and RR flux see Fuji 
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findings, but 17] fails to do so in the sector that involves fermions. However, we also find that 
there are non-zero tree-level amplitudes of the form of two light bosons going to two heavy 
fermions not accounted for in [8| |11| which result in an apparent mismatch at the one-loop level 
(for technical reasons we only check this in the strict a — 1 limit). This can be understood 
in terms of the tree-level amplitudes via the optical theorem, and we explicitly show that the 
additional missing amplitude accounts for the mismatch. This apparent mismatch may be due 
to the a — > 0, 1 limit of [s||TT] being subtle, but from the worldsheet perspective one would 
expect the limit to be smooth. At the one-loop level we also find that the two modified phases 



of [TOJ (one of the phases was first derived in 22 and agrees at least to lowest order) reproduce 



the correct answer while the standard Hernandez-Lopez phase 26 does not, see also [9]. 

The outline of the paper is as follows: In section [2] we review the derivation of the gauge- 
fixed BMN Lagrangian up to quadratic order in fermions [3]. In section [3] we compute tree- level 
amplitudes for two incoming light-bosons. First we generalize the light to light BB — > BB 
processes in |4j to the more general et-gauge, which turns out to be useful in comparing to 
the exact S-matrices of [7 11 . Then we compute additional tree- level processes in the uniform 
gauge (a = |) where we, in particular, find certain non-zero amplitudes not visible in the exact 
S-matrix. In section [4] we utilize the NFS string in the strict a = 1 limit and compute one-loop 
amplitudes for processes with massive and massless in / out-states. In section [5] we compare 
out findings with the proposals of [7j[ll]- We first explain how to match our tree- level results 
with the proposals using some natural choices of the free parameters in [11| together with certain 
gauge-dependent phases. For the one-loop amplitudes we show that the real part of the S-matrix 
element fails to match with either of the proposals for the exact S-matrix in the a — > 1 limit. 
Using the optical theorem we demonstrate that this mismatch is due to the contribution of the 
additional tree-level amplitude not visible in the exact S-matrix. We end the paper with a short 
discussion and outlook and some appendices with details of the calculations. 

2 Action and BMN limit 

2.1 Green-Schwarz action to quadratic order in fermions 

The type IIA Green-Schwarz string action on AdS^ x S 3 x S 3 x S 1 , to quadratic order in fermions, 
takes the form 4F1 



*--,/(!. .% + i . S er A ve - * er.r^e) + o,e% (S .i, 

where the e A (X) (A = 0, 1, • • • ,9) are worldsheet pullbacks of the vielbein one-forms (of the 
purely bosonic part) of the background (* denotes the worldsheet Hodge-dual and we leave the 
wedge product implicit), and O is a 32-component Majorana Spinor. The generalized covariant 
derivative acting on O is given by 

VQ = (V- ^e A T 0129 {l~V)T A ) 9 where V9 = (d - \j ab T ab )Q , (2.2) 

where oj ab is the spin connection of the background space-time. Here V is a projection matrix 
given by 

p = i(i + VSr 012345 + yr^r 012678 ) (2.3) 

and is in fact the projector which singles out the 16 supersymmetries preserved by the back- 
ground. The parameter < a < 1 determines the relative size of the two S 3, s. For unit AdS- 



The IIB case considered in [lj is easily obtained by T-dualizing along the S 1 . 
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radius the S^-radii are 

R+ = -i , R- = -jf= (2.4) 

\/a \J1 — a 

The case a = 0,1 corresponds to AdS$ xS 3 xT 4 where one of the three-spheres is decompactified. 
2.2 Gauge fixed BMIM Lagrangian 

We will use the coordinates described in Appendix A of [4] and we consider a string moving 
along a geodesic involving the angles ip$ and <pg of the two three-spheres, again as in [4]. In order 
to allow for a general light-cone gauge fixing we choose the following light-cone coordinates 

x + = (1 - a)t + a(\/acp 5 + \J\ - a<p 8 ) , x~ = ~(f - (Va<p 5 + \J\ - a<p s )) , 



t = x + + 2ax , v = y/l — a (^5 — \fa (p$ , (2-5) 
f5 = \fa(x + — 2(1 — a)x~) + \J\ — a v , = \/I — — 2(1 — a)x~) — u 

together with 

r+ = T + VaT 5 + y/1 -aT s , F_ = 2(1 - a)F - 2a(^F 5 + y/1 - aT 8 ) , 

where < a < 1 and generalizes the standard uniform light-cone gauge, see [27] for a detailed 
review. The light-cone gauge is then fixed by setting 

x + = r, p + = 1 , (2.6) 

where p + is the conjugate momentum to x~ . The full Lagrangian simplifies considerably in the 
uniform light-cone gauge where a — i and we will mostly restrict to this choice. 



Fixing the gauge and expanding the action (2.1) in transverse fields, with kept fixed, 
gives the BMN Lagrangian 

C = C 2 +g- 1/2 C 3 +g- 1 C 4 + ... 

where g ~ y/X and will be left implicit in most of the calculations. The action describes four 
complex bosons, y\ = ^(^l — ^2) etc, and complex fermions, x± i = 1, . . - , 4. The quadratic 
part is given by [3] (<9± = do ± di) 

£2 = ix\d^ X \ + itf-d+X- + -d+Vid-fc + -d-yid+yt - m 2 yi y % - m^x+x'l + X-X+) > ( 2 -7) 
where the masses are 

mi = 1, m 2 = a, m 3 = 1 — a, m 4 = (2-8) 

and a parameterizes the relative size of the S 3 radii according to ( |2.4| p| 

Thus we have four (complex) coordinates with generally distinct masses. Their charges under 
the various t/(l)'s can be read of in table [l] Note that the massless boson 2/4 has zero charge. 

3 The relation between a and the angle <f> used in |4ll5] is a = cos 2 </>. 
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-1/2 


1/2 


1/2 



Table 1: U(l) charges 



The cubic Lagrangian takes the form 4 



c 3 = 7 =V^ 



l -4 -2 

ot{X-X- 



x-x- + xlxl - t+x 2 +)y2 



- i(l - a)( X 3 _xi + X-X- + X 3 +X+ + X+xV)2/3 (2-9) 

- 2(x-X+ + X+X-)y[ + 2(x-X 2 + - X+X-)2/4 

+ ( x 3 -xX - x 2 -xi)(2/ 3 + y's) + (x 3 +xi - xlx-){h - y' 3 ) 
+ i(f-x\ + xlxDim + y' 2 ) + i{t-xl + xixDih - y' 2 ) 

- a)(a\y 2 \ 2 - (1 - a)\y 3 \ 2 ) y 4 +h.c. , 

where the hermitian conjugate is defined in the standard way, (x-X+)^ = X+X- e tc. Here 
do and d\ are denoted by dots and primes respectively. Note that for the a — and a = 1 
cases, corresponding to AdS 3 x S 3 x T 4 , the cubic Lagrangian vanishes. This simplifies many 
computations in the T 4 case as there are then no 3-vertices. 



In order for the light-cone gauge (2.6) to be valid beyond the cubic approximation, higher 



order corrections to the worldsheet metric are needed. These enter in the form 

1 i 1 2 
7<i = Vij + -7ij + + • • ■ 



(2.10) 



where the subleading pieces are quadratic and quartic in fields respectively. For general a, the 
corrections are rather involved but they simplify in the uniform gauge. For example, the leading 
order corrections for a — \ are 



7oo 



l-|2/i| 2 + « 2 |2/2| 2 + (l-a) 2 |2/3| 2 , 



7f! 



-l-| yi | 2 + a 2 |y 2 | 2 + (l-a) 2 |y3|' 



and = 7 10 = 0. The subleading pieces are involved but straightforwardly determined. Having 
fixed the worldsheet metric so that the equation of motion for x~ is satisfied, one can in principle 
write down the Lagrangian to any order in perturbation theory. Naturally, for general a and a 
the structural complexity is rather daunting. Here we will only write the purely bosonic terms 



4 Note that we are using a slightly asymmetric notation for the 2 and 3 fermions. The action should be 
invariant under 2 «-> 3 together with a «-> 1 — a, a symmetry which is not manifest in our notation. In our 
notation the replacement 2 o 3 should be accompanied by a simple redefinition of the fermions. 
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in the quartic Lagrangian which, in the a = | gauge, take the form 



Cf = a(l - a) (a\y 2 \ 2 - (1 - a)|y 3 | 2 ) 2 - ~a(l - a) $ + # - - (|y 2 | 2 + |y 3 | 2 ) 



- |y 4 | 2 (| yi | 2 - (2a - l) 2 (a|y 2 | 2 - (1 - a)|y 3 | 2 )) + \yi\ 2 (a 2 \y 2 \ 2 + (1 - «) 2 |y 3 | 2 ) 

- (I2/2I 2 + |y 3 | 2 + Ml 2 ) (M 2 - a 2 \y 2 \ 2 (1 - a) 2 |y 3 | 2 ) - afe| 2 M 2 - (1 - a)N 2 |y 3 | 2 
+ a(l - aM| 2 (M 2 + I2/3I 2 ) - lyiHyil 2 + «|y 2 | 2 |y 2 | 2 + (1 - «)|2/ 3 I 2 M 2 . 

(2.11) 

The relevant terms for our computations which mixes bosons and fermions can be found in 



(E.l) 



3 Tree-level scattering in AdS 3 x S 3 x S 3 x S 1 



Following the notation of 28 , the worldsheet S-matrix separates into scattering, transmission 
and reflection pieces, 



S = 1 + -S- 
9 



T = 1 + -T 
9 



R = -R 
9 



= 1 
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(*) 



2 ( S ) 3 



X 

(c) 

X 

(c) 

X 

(*) 



We will consider processes with two incoming bosons and we write the relevant S-matrix elements 
as 



T • \yi{pi)yM) = A^ k % k y{) + B^^\ X kXi) 

R • \vi(pi)yj{P2)) = A (l3M) \ym) + B {i ~ m \xkXi) 



(3.1) 



where we will take i,j = 2,3 to be any of the two light particles with masses a and (1 — a) 
respectively. When massless particles are involved a somewhat more general structure is allowed 
with A l i' kl etc. Note that many amplitudes are trivially zero due to [/(l)-invariance. 

It will turn out that for all the amplitudes we compute the i?-piece will be zero and we will 
not write it out explicitly. 



Diagonal light processes in general a-gauge 

In [i] the diagonal elements — for i, j — 2,3 were determined in the uniform light- 

cone gauge a = h. We will begin by generalizing this result to a general a-gauge. Summing up 
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the s, t and it-channel contributions together with the contact terms we find 



= 1 + \ ^ +P i + Ul- 2a)(J* P2 - J? P1 ), (3.2) 

Z LU\ 'p 2 - 'Pi Z 

4(22) , , i a(pi-p 2 ) 2 i , (2 ) (2) \ 

A{ ' = 1 + ^-7^ m~ + -I 1 - 2a) (wi 'p 2 - w 2 Pi) 

A uj[ 'p 2 - w 2 'pi L 

A (2 3) = A m = l + 1 (1 _ 2a){^p 2 - 4% • 

The 33, 33, 23 and 23 elements are given by sending a -> 1 - a in an obvious way. We have 
used the fact that the delta functions from energy momentum conservation forc^] 

P3=Pl, PA=P2 (3.3) 



and the relativistic on-shell energies above are explicitly given by 

wW(p) = v / l+p 5 , uj (2) {p) = V" 2 +P 2 , ^ (3) Cp) - Va-c^+F 2 , ^ 4 )( P ) = |pI 

In the above equations we have included the external leg factors and Jacobian from the overall 
delta functions, 



A(Lu[ l) p 2 ~4 l) p 1 )' 

We now turn to the non-diagonal elements in the scattering of two light bosons, as some of these 
are rather involved we will only evaluate them in the uniform gauge a = |. 

Non-diagonal processes in uniform gauge 

Here we will consider processes where the incoming particles are light-bosons but the out-going 
particles can be anything allowed by conservation of the {7(1) charges. For transition type 
amplitudes, where the in-state is neutral, there are typically several possible out-states. What 
is more, when the masses of the scattered particles are different the solution of the energy 
momentum constraint becomes more complicated. Since some amplitudes are rather involved 
already at tree-level we will restrict to a = 1 in these cases. Collecting the various contributions 
we havfS 

\V2V2) -> A^\ Vl y ) + B^\ Xl Xj) : 



^(22,11) = q ^(22,33) = q ^(22,44) = ^(22,44) = ^(22,44) _ q _ 

5(22,11), =i = 1 (1-P 2 -)(1-(P'-) 2 ) ^ B( 22,22) = i (a 2 -pl)(a 2 -(p>- f) ^ 
4 y/p_jj_(p_ +p'_) ' 4a y/p_p'_(p_ +p'_) 

5(22,33) _ ^(22,44) = ^(22,44) _ £(22,44) = q ^ /ggN 



5 With particles of different masses, this relation can become significantly more complicated. 

6 To see that processes involving particles of different mass, for example A^ 22,33 ' and A^ 22 ' 11 ), are zero a 



useful identity which follows from energy-momentum conservation is 

m 2 (wi^ 2 + P1P2) = m 2 (w 3 w 4 +P3P4) , 
where the incoming particles have mass m; and the outgoing ones mass mf. 
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Figure 1: The s,t and it-channel diagrams. We also have six-vertex contact diagrams but these have 
trivial topology 

where p_ = p\_ and p'_ = p 2 _ and the light-cone momenta satisfy p + — ^- since the particles are 
on-shell. Later we will utilize the so-called Near Flat Space, or Maldacena-Swanson, limit 24 
where the right moving sector is boosted and we choose the right moving notation here to 
establish a simple connection with the results presented ther^] In the mixed sector we find 



#(23,23) 



{a 2 ~pl){{l- a f-{p'_f) 
A^Jp_p'_ ((1 - a)p_ + ap'_) 



B (23 ' 14) Ui 



lp'-(l-p 2 -) 
4 ^/p—fT 



(3.6) 



Let us also note that in the purely bosonic sector there is no mixing between the light and the 
massless modes at tree-level, i.e. 



^(22,44) _ ^(22,44) = ^(24,24) = ^(24,24) = q ^ 



(3.7) 



This follows, for example, from (3.5| by crossing 



4 One-loop scattering in AdS s x S 3 x T 4 

Generally we have both cubic, quartic and sixth order vertices implying that the topology of a 
generic one- loop amplitude can be rather involved. In order to reduce the number of contributing 
topologies we will work in AdS 3 x S 3 x T 4 , i.e. a = 1, and take the a — | gauge . In this limit 
the third (complex) coordinate and x± become massless and the cubic Lagrangian vanishes. 
Equivalently we could of course take a — where ?/2 and x± become massless. In both cases 
the contributing diagrams are just of four distinct types, see figure [T] 

4.1 One-loop amplitudes in the massless sector 

In this section we will compute the S, T and R matrices for massless bosonic 2 — > 2 processes. 
For these modes the dispersion relation is simply e = \p\ and as it turns out the final answer 
will take a rather simple form. Furthermore, the relevant part of the quartic BMN Lagrangian 
simplifies to 

U - (4.1) 

-i (d + y 3 d_y 3 + a + y 4 a_y 4 ) {t+xt + xlxl) - \ (d-y 3 d + y 3 + d + y 4 d„y 4 ) {xtx% + t + X-) 

-~(d + y 3 d + y 3 + d-y 3 d-y 3 + d + y 4 d + y 4 + d_2/ 4 <9_?/4)(M 2 - I2/2H 

-^(d^y 4 d + y 3 - d+y 4 d-y 3 ) (x+xi - - l,(d+y 4 d~y 3 - d+y 3 d^y 4 ) (x%X- ~ xtx\) 

7 A11 expressions presented in this section are however for the BMN string, albeit in a right moving notation. 
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Thus the one-loop diagrams with all massive (or massless) external particles will have either two 
massive modes or two massless modes propagating in the loop. 

The full sixth order BMN Lagrangian is naturally rather involved, but as it turns out we 
only need the purely bosonic part for massless processes, 



Aj = -|(<9+?/3<9 + 2/3 + d + y 4 d+y 4 ) (d-y 3 d-y 3 + d-y 4 d^y 4 ) (\ 



l/i I 



12/2 1 



(4.2) 



Furthemore, since the particles are massless it makes sense to demand that the sign of the 
incoming momenta should be opposite, otherwise they would never meet. Hence we will assume 



Pl > > p 2 . Evaluating the integrals in ( C.l )-(C.3 ) using (B.l ) with the notation p± — p± and 



P± = P± S 1VCS 



? (i) . 

'33 

i 



T (l) 
33 

i 

+ 8n 



i 

8^ 



i 

8^ 



{p+p'-) 
- - + 7 + log7r)(p + p'_) 



2 \ 2 

- + 7 + log 7T - \0gp'_p + J (p+p'_) 

2i 

= — (l-log [- 4piP2])Pipj, 



I 

8^ 



[P+P- 



i 

87 



2 

-+7- 



■ log 7T - log -p'_p+ ) (p+p' 



(4.3) 



(4.4) 



+ 7 + log7r)(p+p'_) 



2i 



{l-\og[^p 1 p 2 ])p\pl, R { £ = 0, 



Sil ] = * (l - log [A PiP2 ])pIpI T$ = |(l - log [ - A PiP2 })pIpI 



? (1) 



2i 



~ (1 - log [A PlP2 ] )p\p\, T 4 ( 3 X) = ^ (1 - log [ - A PlP2 ] )p\pl, = (4.6) 



2i 



R<V = 0, (4.5) 



where the c index denotes contact terms arising from six-vertex interactions. We have ignored 
the external leg factors and Jacobian from the overall delta functions. Furthermore, note that 
in the NFS limit, see (A.2) for details, the amplitudes are zero since p+p'_ ~ 0(1). 



Since P\p 2 < (or equivalently p_ < p'_ ) we see that some of the processes develop a real 
part which breaks the 5*0(4) invariance. This should be due to the fact that we are working 
in type IIA where there is flux along one of the T 4 -directions. T-dualizing the x 9 -coordinate 
one should recover the (manifest) 5'0(4)-invariance of the type IIB case. We will leave further 
investigation of the massless processes for an upcoming paper. Especially interesting is of course 
to understand the apparent mismatch with the conjectured set of Bethe equations for massless 
fields in JfjJ. 



4.2 Massive one-loop scattering in NFS limit 

We now wish to look at the scattering of the massive bosonic modes in AdS 3 xS 3 xT 4 (a = l). 
However, at the one loop level one encounters problems with the BMN string since the am- 
plitudes (naively) fail to be UV finite, a problem which also appears for the more symmetric 
AdS$ x S 5 string This is most likely related to some subtlety with the light-cone gauge fix- 
ing we employ. Luckily, the problem can be circumvented by utilizing the Near Flat Space, or 
Maldacena-Swanson, limit [24]. The details of the limit can be found in appendix A, see espe- 



cially (A.2). The limit basically boils down to a boost of the right-moving (worldsheet) sector 



and all divergent loop integrals are shifted to subleading orders in the l/g expansion. Further- 
more, simple power counting demonstrates that only quartic vertex diagrams contribute, so the 



8 We thank T. McLoughlin for making us aware of this (unpublished) result. The reason this problem did 
not appear already in section |4.1| is probably because there only the purely bosonic piece of the sixth order 
Lagrangian contributed. 
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topologies of Feynman diagrams are completely summarized in figure [T] 

In a = 1, the massive modes are the 1 and 2 particles. Using the expressions in (D.l) we 
find0 



i (p-p'-Y 

32tt {p'_) 2 -p 







/ \2" 




( 


i(p_ 


+ P-) 


+ 

t L 



i(p'_ + P-)((p'_) 2 +p 2 _) . pL] r(p'_+P-) 



pi - p_ 



log 



P- 



(7r + i log — 

pL — P- P-' 



i (p-p'-Y 

32tt (p'_Y - p 2 _ 
(i) _ 1 (P-P- 



32tt (p'_) 2 -p 2 . 

2 



i (p_ + p' 

2 

X 



, 2 2ip_p'_ (P- + PL ) log ^ tt ( p _ + p'_ ) ; 



pL -p_ 



p'_ -p_ 



(- «(p'--P-) 
1 (p-p'_) 2 



j(p'_-p-) p'_ 
, T ~ l og — 

p'_ + p_ p_ 



r ( p , --P-)((p , -) 2 +P 2 -) / p'_. 

j^- '- (tt + i log — 

pL + P- P- ' 



(i) 



32tt (p'_) 2 - pi 

1 (P-P'_) 2 
327r(p'_) 2 -pi 

I \2 



HP- -P-; 



HP--P-J + 



2ip_p'_ (pi - p_ ) log ^ 7r(pL - p_ ) (p± + (pL ) 2 ) 
pL + p- p- + p'_ 

2ip_p'_(p'_ -p_)log^ 



■ p- 



T (i) = J- (P-P-) 
y 327r(p'_) 2 -p 2 _ 



HP-+P-) - 



+ 7r((pL) 2 -p 2 _) , i±j 
p_ + p_ J 

2ip_p'_ (pi + p_ ) log ^ 7r( p _ + y_ ) (pi + (p'_ ) 2 ) " 



PL - P- 



pL - p_ 



(4.7) 



i 7^ i 



where we, in contrast to the massless case, also included the external leg factors and the Jacobian 
from the delta function, which in right moving components equals 

1 p_p'_ 



2{p'_) 2 -p 2 _ ' 

While we have not written out the explicit contributions, one can also show that the s,t and 
u-channel contributions in the reflection piece, cancel among themselves. Hence, at least 

for a — 1 the one-loop S-matrix seems to be reflectionless, a fact so far only seen at tree-level 
in (4). 



5 Comparison to proposed exact S-matrix 

Here we will compare our results from worldsheet calculations to the proposals for the exact 
S-matrix 117]. Since we are looking only at the scattering of bosons the relevant terms in the 
S-matrix are 

S\ Vi v j )=AM\y i v J )+BM\xiXj). (5-1) 



3 We have taken out a factor of l/(167r) and (p'_P— ) from the s,t and u-channel contributions in the first 



lines of and t£\ 
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For the S-matrix proposed in [TT] these take the form 



together with 



B 



(22) 
12 



B 



(23) 
12 



(22) _ 



A 12 
(22) 



12 



4(23) 
A 12 



A 



(23) 
12 



(22) n 



rc (22)N-l^l - ^2 

l°i2 ; ^+ _ 



(22)x-l 



(sf 2 3) )- x 



(i-^) (i-^ 



x iv2 



_j -^l")^ -^2") 



_j y-(^i -^i")^ -2/2") 



^1 y 2 



1 



x tvi 



1 



(5.2) 



(5.3) 



where particle 2 and 3 have mass a and 1 — a respectively. The other elements are obtained by 
2 o 3 and x O y. The phases appearing in these expressions are given by 



4,(22) _ ^(22) 
D 12 — °12 



1 



4,(23) _ £,(23) I z+s/ 2 
°12 — °12 



1 - 



x i vi 



1/2 



= r -^(^iP2-" 2Pl ) x l x 2 



1 



3v ^ X 2 1 



■^2 2/ \ 



-1/2 



x{ y 2 / 



1 



x i vi 



(5.4) 



Note that we have chosen the free parameters 7 and T appearing in [TT] to be zero so as to 
have the standard factor of a 2 in S^ 22 ) that one expects to have [lj. In [jj S^ 23 ) was taken to 
be simply 1, corresponding to T — —1/2 in [TTJ , but this does not seem to give the right answer 
in a general a-gauge as we will discuss below. We have also translated to the string basis by 

f x - x+X 1 / 2 

multiplying with I -?r ) . For this choice of parameters the purely bosonic S-matrix also 

V x _\ x 2 / 

agrees with that of [71, modulo the issue with S^ 23 \ provided that one makes the replacement, 



x t x i) x i x 2 ) 



\ 



1 - 



(5.5) 



in A^\^ and similarly in For the S-matrix of [7] the B coefficients vanish for a 7^ 0,1. 

This clearly disagrees with the results of our worldsheet computations since wc find B 7^ 0. 



Note that we have included in (5.4 1 gauge-dependent phases [27| involving parameters a and 



b, which will be related to the worldsheet a-gauge parameter below. These gauge-dependent 



phases were not included in the analysis of 11 but we seem to find that they are needed in 



10 



order to match our worldsheet results in a general a-gauge. 

The phase a appearing in (5.4) is required to reduce to the (modified) AFS phase, discussed 
below, at lowest order. 



5.1 Tree-level 



Here we will expand the S-matrix elements in (5.2) and (5.3) to lowest order in h(X) 



+ 



0(A ) and compare them to the worldsheet tree-level results found earlier. Using the fact that 



^Ja 2 + 4/i(A) 2 sin 2 J 



2h(X) sin -f 



2A ±i 

e 



2 A 



1 ± i 



(5.6) 



and similarly for y ± with a 



A 



(22) 
12 



A (22) 
A 12 



= 1 



1 — a, we get 

i a{pi+P2) 2 
2h lo\P2 - 

i a(pi-p 2 ) 2 
2h LO1P2 - U2P1 



2ha 

i 

2ha 



(1 - 2a)(u>ip2 - LO2P1) 
(1 - 2a)(u>iP2 - UJ2P1) 



A 



(23) _ ,4(23) _ 



A)i' = 1 - 



.1-25 
2ha 



-UJ1P2 + i 



1 - 2b 
2/i(l - a 



-UJ2P1 



(5.7) 



This agrees completely with the tree- level S-matrix computed in the near BMN- limit (3.2) 
provided that a and b are related to the a-gauge parameter as 



S=i(l-a(l-2a)), 



6=i(l-(l-a)(l-2o)). 



(5.8) 



Unless a = \ these cannot both be set to zero by a choice of a. 

Note that for the S-matrix proposed in [7] S^ 23 - 1 = 1 which would give A^ 2 ^ = 1, disagreeing 
with the worldsheet calculation except when a = | . Modulo this issue and the gauge-dependent 
phases we see that the S-matrices of [7] and [Tl] match the worldsheet calculation in the bosonic 
sector at tree-level. 



For the scattering of bosons into fermions we get for the S-matrix of 11 by expanding (5.3 1 
(setting a = \) 



D 



(22) 

12 



wi)(a - lo 2 ) 1 



D 



(23) 



2h 



\J(a - wi)(l - a - L0 2 ) ( 1 - 



Pi -P2 
W1P2 - W2P1 

(1 - a)pi 



1 (a 2 - p 2 _)(a 2 - (p'_) 2 ) 



4ah 



ap 2 



U1P2 - U2P1 



\/P-P'- (P- +P-) 
1 (a 2 -p 2 _)((l-a) 2 -(pi_) 2 ) 
4:h ^p_p'_ ((1 - a)p_ + apL) 

(5.9) 



These expressions match exactly (apart from a factor of i which is due to our conventions for 



the fermions) those coming from the worldsheet calculation in (3.5) and (3.6). 

In conclusion we find complete agreement at tree-level between the worldsheet computations 
and the results coming from the proposed exact S-matrix of 11 with gauge-dependent phases 



added. For the amplitudes involving fermions we disagree with the S-matrix of [7]. 



5.2 One-loop comparison 



To lowest order the phase a appearing in the S-matrix is the AFS phase 29 . Its standard 



all-order extension, the BES phase, does not seem to give the right answer at the one-loop 
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order however [30]|9] . Instead a proposal for the one- loop correction to the phase was presented 



first in [22] and then improved in 10 . Here we will expand this phase in h(X) in order to 
compare to the one-loop NFS computation on the worldsheet. Since this computation was done 
in AdS^ x S 3 x T 4 only we will take a = 1 here as well. The full phase has an expansion^] 



a(k,j) = exp t6kj 



1 kj 



1 

2m 



E E 



r— 1 s— r+l,r+s odd 



qr{xk)q s {xj) - q r {xj)q s {x k ) (5.10) 



where m — (a, 1 — a) for particle 2 and 3 respectively (in the mixed sector the phase takes a 
slightly different form) ancj^j] 



1 



1 



71- 1 \{x+) n ~ 1 (x-)™- 1 

Expanding the magnon charge 



c r . s = h(X)S r+1 . s + —cl 1 } + 
zir 



( , Pk (u k - ax"- 1 _ 3 

q n {x k ) = -r[ +0(h °) 

h \ p k ) 



(5.11) 



gives to leading order the AFS phase 

h 



nAFS 



{xk j Xj ) 



2a 



y^[gr(fffc)gr+l(ffj) - <7r(£j)gr+lOEfc)] 
r=2 

1 

4ft 



\ , (Pj - Pk) . 1 , V 

2(Fj - Pfc) + a h - {ujjPk - WfePj) 



+ 0(fr- 3 )(5.12) 



For the mixed (23)-sector we use instead the form of the phase given in [XT] 



qAFS 



(^fe ; %) = log 




1 

2h 



(Pj - Pk) 



8ha(l — a) 



1 (1 - a)p k - apj 
((1 - a)p fc - apj) 2 

PjUJ k - p k Wj 



+ PkUJj - PjLOk 



+ 0{h~ 6 ) 



(5.13) 



Since the AFS phase is exact to order 0(h~ 3 ) the one-loop contribution comes only from ci]} 
together with sel f int eractions of the tree-level contribution. In the NFS limit we find for the 
phase defined in (5.4) that only the one-loop piece of a involving ci]l contributes 



-,(22) 
'12 



(5.14) 



Since this gives only an imaginary contribution the real part of the one-loop amplitude A^ 22 ^ 
should therefore, according to 11,8] eq. (5.2), be given only by 



(p-P'-fiP- +p'-f 



32 (p_ -p'_ 



+ o(h- 4 ) 



(5.15) 



10 The sum starts from r = 1 instead of r = 2 to account for the form of the one-loop phase of [10| 
xl The factor of ~ in front of the one-loop term follows from our definition of the coupling. 
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where the ellipsis denote tree-level terms. This is in complete agreement with the real part found 



in (4.7). 



Similarly we find that the real part of A^ 22 "> in the NFS limit comes from 



1 - 



1 



1/2 



for the S-matrix of 11 and 



P-P- 



{p'-f 



32(p_+p'_) : 



(5.16) 



(p-p'_) 2 (p--p'-¥ 

32 (p_ +p'_) 2 



(5.17) 



for the S-matrix of 7 , neither of which agree with the real part of in (Wjl. The reason for 
this disagreement will be discussed further in the next section by applying the optical theorem. 
Let us now look at the imaginary part of the amplitudes. As mentioned these can only receive 



contributions from the one-loop phase. In 10 the form of the one-loop contribution to the phase 
was suggested to be 



s + r - 2 ' 



s + r - 2 



(5.18) 



where bar indicates T-type scatterings. Switching from the k and j labels to p and p' a somewhat 
lengthy calculation gives the NFS expansion^] 



iO(p,p') 



/\ 1— loop 



1 



2n r + s - 2 



q r {x)q s {x') - q r (x')q s (x) 



(5.19) 



i0(p,p') 



f\l— loop 



% (P-P- 
64tt (p'_Y - 

2tt ^ 



P- 



2p-p'_(p- log 

p'_ - p_ 



p- 



P- 



o(h- A ) 



s — r 

\ 2 



q r (x)q s (x') - q r (x')q s (x) 



[P-P- 



64tt (p'_) 2 -p 



2p-P-(p'- -P-) log 



V- 



P- 



(P+ -P-j 



0{h- 4 ) 



The imaginary ter ms i n the amplitude are (—2) times these expressions since the phase enters 



as a 



and Tjp in (1471) 



in (5.2). We see that these expressions nicely match the imaginary parts of S, t 



(i) 



Thus we find that the expressions coming from the S-matrix proposed in [Til for 22-scattering 



at a — 1 and in the NFS limit match with those coming from the worldsheet calculation in (4.7) 



except for the real part of A^ 22 \ Furthermore, to get the matching of the imaginary parts we 
needed to use the one-loop phases proposed in IIO) instead of the standard Hernandez-Lopez 



phase 26 . We will now analyze further the reason for the discrepancy in the real part of A 1 - 22 ' . 



5.3 Optical theorem 



Since the disagreement we found in the last section with the proposed S-matrix of 11 



is m 



the real part of the S-matrix element (imaginary part of the amplitude) we can use the optical 
theorem to trace it back to the tree-level. Taking a = 1 and the NFS limit in the tree-level 



2 The first of these agrees with the phase obtained in 22 
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results in (3.2) and (3.5) we find 



\V2y2) H- (1 + A^ 2 y 2 )+B^\ XlXl ) + B^\ X2 x 2 ) : 

(5.20) 

A /(22) = _ • P-P-(P- -P-) B (22,ll) = ^(22,22) _ _ * (P-P-) 3/2 

4 P- + 4 p_+ p'_ 

The optical theorem then tells us that the real part of the S-matrix element for 22-scattering 
should be given by 

|^(22)|2 + | B (2241)|2 + | B (22,22)|2 = Up-T>-?(j£ + 

in agreement with what we found for T 22 in (|4.7|) (the relative factor of (—2) is due to our 



conventions). If we would omit the amplitude for y 2 y 2 — > XiXi we S e t instead the result of 11 



in (5.16) 



The process y 2 y 2 — > XiXi is n °t accounted for in the S-matrix of [IT] because for a ^= 0,1 the 
fermion xi is heavy and therefore treated as a composite object (28,5 . Our calculations show 
however that there is a non-zero amplitude for the process y 2 y 2 — > XiXi at tree-level and this 
affects other one-loop S-matrix elements through the optical theorem. It is clear then that the 
proposed S-matrix of [IT] should somehow be modified in order to take this fact into account, 
at least in the case a — 0, 1. Indeed at a = 0, 1 the j/i and xi excitations are no longer 'heavy' 
but become fundamental from the Bethe ansatz perspective so one expects the exact S-matrix 
to include also these processes. This should resolve the one-loop discrepancy at a = 1, but from 
the worldsheet perspective one expects the S-matrix to change continuously with a so it is hard 
to believe that the discrepancy would not persist for a ^ 0, 1. 

6 Discussion and outlook 

We have computed tree-level worldsheet S-matrix elements for the string in Ad S3 x S 3 x S 3 x S 1 
and certain one- loop elements in AdS% x S 3 x T 4 . The tree-level results are in agreement with 



those coming from the exact S-matrix proposed in 11 provided one includes appropriate gauge- 
dependent phases. The tree-level amplitudes involving fermions disagree with those predicted 
from the S-matrix of (7j. 

At the one-loop level we find agreement for the imaginary part of the S-matrix elements we 
computed provided one uses the one-loop phases of (lO] (one of which appeared in [22]) which 
differ from the standard Hernandez-Lopez one [26] . This gives a nice check of their results which 
were derived by different methods. For the real part of the S-matrix element we seem to find a 
disagreement with the proposed exact S-matrices however. We have traced this issue to the fact 
that at tree-level there is a non-zero amplitude for two light scalars going to two heavy fermions 
which is not visible in the exact S-matrix, since for general a the heavy modes are thought of as 
composite. This tree-level amplitude is related to the mismatch in the real part of the S-matrix 
element at one-loop level through the optical theorem. This mismatch could be due to a subtlety 
in taking the a — > 0, 1 limit in the S-matrix of [11| . Indeed one expects that at a — 0, 1 the modes 
that are heavy for < a < 1 should no longer be composites but now appear as fundamental 
excitations and the S-matrix should take this into account. It would be very interesting to see 
whether there is also a one-loop mismatch for a ^ 0, 1 but this would require considerably more 
complicated calculations. From the worldsheet point of view one would expect this to be the 
case since the limit a — > 1 is completely smooth, at least at the level of the Lagrangian. 
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We have also computed several S-matrix elements involving the massless fields. It would also 
be interesting to understand better how these are accounted for by the Bethe ansatz. We hope 
to return to this question in the near future. 

Finally it would be interesting to carry out similar calculations in the less supersymmetric, 



and much less studied, case of Ad5 2 /Ci 7 'T 1 [31 33 
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A Details of the near flat space limit 

The near flat space (NFS) limit or Maldacena-Swanson limit 124] can be considered as a world- 
sheet Lorentz / boost transformation together with a resummation. Originally the BMN La- 
grangian has the form 

Cbmn = £2 + -^£3 + + (A.l) 

where the subscript denotes number of transverse fields in each term. Performing the following 
worldsheet transformatiorF^l 

p±^g T h±, d±-> g T ^d±, x ±->g Tli x± (A.2) 

in the BMN Lagrangian and taking the g — > 00 limit gives a theory without coupling dependence 
but still exhibiting non-trivial interactions 

£ N FS = £2+£3 + £4 + 0(g~ 1 / 4 ) (A.3) 

where the higher order terms £3 and £4 are considerably simpler than the corresponding ones in 



(A.l ). Furthermore, all vertices beyond quartic order are suppressed in the NFS limit, simplifying 



it further compared to its BMN counterpart. 

B Integral conventions 

The integrals we encounter are standard ones 



B„ 1 ...„ p [pt,P2,mi,m 2 \= I -iXi-T^ ^fr"^^, ^, (B.l) 



(27r) 2 (fe2 _ m 2) _p L _ p^2 _ m 2) ' 

T \m]= [ k ^--- k ^ 
J (2 ^)2 ( £ 2 „ to2) 

where the B and T label denote bubble and tadpole type diagrams. 



13 In order to make contact with the scalings of [24] one should remember that the BMN has already scaled 
the fermions with a factor of g~ 1 ' 2 . 
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NFS integrals 



The integrals we encounter for the NFS string are versions of (B.ll with pure powers of fc_ in 
the numerator. Furthermore, for processes with massive external legs the two particles in the 
loop always have the same mass. Using this and denoting the power of fc_ with r we can make 
a slight improvement of the notation as 



B r [p,p',m] = J d 2 fc- 



(fc 2 — m 2 ) ((fc — p — p') 2 — m 2 ) 
which after introduction of a Feynman parameter becomes 



B r [p,p',m] = ^ J 



dx ((l-x)(p-+p'_)) r 
m 2 + (x — l)x(p + p'j 



Assuming that p\ > p[, which implies p'_ > p_, and using an IR regulator p we find for the 
relevant integrals 



s-channel: 



1 P-P- P- , ^ n r- -/ ,1 1 P-P 1 - P- , \ 

7T 



B [p,p',l] =-^-TTT2 r(* lo S h7r )> B i[i>,p'A] = -7--, ( ilo & 

1 J ZTT [p_Y —p_ P- 47T p_ — P- jj 

u- channel: 

Si [p, -p', 0] = ±- J~ P p _ log P, B 2 [p, -f, 0] = - ^P-P- (1 + log p) , 
B 3 [i?,-p',0] =^p-p'_(p'_-p-)(Z + 2logp), S [p, -f, 1] = ± J_yZ p 2_ Jog ^ 

i p _ p >_ P '_ ^ ^ i P -p'_((i-iog^)(y_) 2 -(i + iog 



Si [p, -p", 1] = log^, S 2 [p, -p', 1] = - 



Air p'_ + p- p-' 4ir (p'_) 2 — P- 

We also have t-channel processes. These however turn out to be rather subtle when we have 
massless particles propagating in the loops. The best way to evaluate them is to perform the x 
integration using the u-channel integrals and at the end expand around p' — > p. 
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C Details for massless scattering 

The one-loop elements for the massless processes sQ , and R^' , are encoded in the following 
expression}^] 

Sg } = 2(e 2 +p 2 )(e 2 +p 2 )B[f 1 ,-p 1 ,l] - (p 1 ) 2 {p 2 ) 2 B + _[ft,-ft,0] (C.l) 
+ (2(ft • ft) 2 - (ft) 2 (ft) 2 ) ((A - ft)+B_ [ft, -ft, 0] - B + _ [ft, -ft, 0] ) 
+2( £l e 2 + PlP2 ) 2 B[p u -ft, l] + 2(ft • ft) 2 T[l] , 

= 2(ef +p 2 )(e|+p 2 )B[ft,-ft,l] -(pi)2(^)2s + _[ft,-ft,o] 
+ (2(ft -ft) 2 - (ft) 2 (ft) 2 )((ft + ft)+B_[ft,ft,0] 
+2( £l e 2 +pip 2 ) 2 J B[p 1 ,p 2 , 1] + 2(ft • ft) 2 T[l] , 

=2(e 1 e 2 + PlP2 ) 2 B[p 1 ,p 2l l] + (ft ) 2 (ft) 2 ((ft + ft) + B_ [ft ,ft, Cl] - B+_ [ft,ft, 0] ) 

+(ft) 2 (ft) 2 ((ft - ft) + -B- [A. -ft. 0] - B+- [ft, -ft, 0]) 
+2(e 1 e 2 + :Pl p 2 ) 2 B[ft,-ft,l] + 2(ft) 2 (ft) 2 T[l] 

while the scattering among the second set of massless coordinates gives 

=2(e 2 +p 2 )(e 2 2 +p 2 2 )B[p 1 ,-p 1 ,l]+2(e 1 e 2 +p 1 p 2 ) 2 B[p 1 ,-p 2 ,l] (C.2) 
+2(ft -ft) 2 ((ft+ft) + B_[ft,ft,0] -B + _[ft,ft,0]) + 2(ft-ft) 2 T[l], 
T$ =2(e 2 +p 2 )(e 2 +p 2 )[p 1 ,~p 1 ,l] + 2(e,e 2 + Pl p 2 ) 2 B[p u p 2 , l] 
+2(ft-ft) 2 ((ft-ft) + 5-[ft,-ft,0] -B + _[A,-ft,0]) + 2(ft .ft) 2 T[l], 

fl£> = -2(ft) 2 (ft) 2 £+4ft,-ft,0] +2(e l£2 +p 1 p 2 ) 2 (i?[p l! p 2! l] + B[ft,-ft,l]) 
+2(ft) 2 (ft) 2 T[l] 

For processes mixing different particles we find 

S$ = 2(e 2 +p 2 )(e 2 +p 2 )B[p 1 ,-p 2 ,l] + (2 (ft -ft) 2 - (ft) 2 (ft) 2 )T[l] (C.3) 
+ 2^1 -Pi) 2 (^2+p 2 ) 2 ((pi -ft) +J B_[ft,-ft,0] -B+_[ft,-ft,0]) 

+ 2^1 +Pi) 2 (e2-p 2 ) 2 ((pi + ft) +J B_[ft,ft,0] -.B + _[ft,ft,0]), 
T$ = 2(e 2 +p 2 )(e 2 + p 2 ) 5^,-^2,1] + (2(ft -ft) 2 - (ft) 2 (ft) 2 )T[l] 
+ 2^1 ~Pi) 2 (e2+p 2 ) 2 ((pi + ft) + 5-[ft,ft,0] -B+_[ft,ft,0]) 

+ o( £l +Pl) 2 ( £ 2 ~-P2) 2 ((pi -ft) + 5_[ft,-ft,0] --B + _[ft,-ft,Q]), 

i?^ =2(e 1 e 2 +p 1 p 2 ) 2 [p 1 ,p 2 ,l] - i(ft) 2 (ft) 2 B+-[ft,-ft,0] 

+ ^(ft) 2 (ft) 2 ((ft - ft) + S_ [ft, -ft, 0] - B+_ [ft , -ft, 0]) + (ft) 2 (ft) 2 T[l] 

14 Expressions like (pi ±P2)_ etc mean the minus component of the vector p\ ip2- 
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D Details for massive scattering 



Focusing on y 2 particles and using (B.ll we have, 

S22 = \ [p-P- (p-p'-Bo [p, -P, 1] + B 2 [p, -p, 1] - B 2 [p, -p, 0])] t (D.l) 
+ \[(p^+p'-) 3 ((p-+p'-)B [p,p'A]~B 1 [p,p',1]^ s 

+ \ [(pi + P 2 - (P-f - P- (P-f + (P'_ ) 4 ) Bo [p, -p 1] + (p- - pi ) (2p 2 _ - p_p'_ + (p'_ ) 2 ) Si [p, -p' , 0] 
-2(p 3 _ + p_ (p'_) 2 - (pi) 3 )Bi [p, -p', 1] - (4pL - 3p_p'- + (p'_) 2 )B 2 [p, -p', 0] 
+ (pl + {p'_f)B 2 [p 1 -p\l] 



T (i) _ _ 
^22 — 



(pl + {p'_f)B 2 [p,p',l] 



p- - pi) 3 (Si [p, -p', 1] + (p'_ - p-)B Q [p, -p', 1] ^ 
-2(p 3 _ + p_(pi) 2 + (p'_f)B 1 [p, I /, 1] + (/_ + p 2 _(pi) 2 +p-{p'-f + (pi) 4 )Bo[p,pil] 
+2p_B 3 [p,p",0] - (4p 2 _ + 3p_p'- + (p'_) 2 )B 2 [p,it',6\ + (p_ +p'_)(2p 2 _ + p_p'- + (p'_) 2 )i?i[p*,p",0] 

+1 [p_p'_ (B a [p, -p, 0] + £ 2 [p", -f, 0] ) + (p_p'_) 2 (B [p, -p, l] + B [?, -f, l] ) 

where we made use of the trivial fact that t-channel integrals with r > 3 are zero and B r [p, —p. l] 
for r = 1, 2 also equals zero. For the remaining 33 and 23 processes we have similar expressions. 



E Relevant piece of quartic Lagrangian 

Here we collect the terms in the quartic Lagrangian of the form (yf , y 2 y 3 , y 3 ) x x 2 

[(xixi + (xi)'xl) - 4i (1 - a) (xixi - xix^)] I2/2I 2 

r(l - «) 2 [(Xixi + (Xi)'xi) + 4ia (xixi - xixi)] 1 2/3 1 2 



/*4 '2 



+ 2 ( a xixi + X 2 +X- + (1 - «)xtxi) 2)22/2 - -£*X-X- 2/2(2/2 - aj/a) 

i 1 

- ^axiXV 2/2(2/3 - aij 2 ) + - ((1 - a)xVxi + xix+ + a***!) 2/32/s 

- i(i - a) (xixi - xixi - xixi + xixi) 2/3(273 - (1 - a) 

- 1(1 - a) (xixi - Xixi - Xixi + Xixt) 2/3(2/3 - (1 - a) 27 3 ) 

- ixixi 2/22/3 - Ixixi 2/ 2 2/3 + 7,xixi ((1 - a) 2 2/32/2 - a 2 y 3 y 2 ) 

- ^Xixi ((1 - a) 2 y 3 y' 2 - a 2 y' 3 y 2 ) + \{at-t+ ~ (1 - «)xixi)2/32/ 2 



+ l(«xixi - (1 - «)xix 4 +)2/22/;3 + \{o?y' 3 y 2 (1 - a) 2 2/ 2 2/ 3 )xi^ 
+ ^((1 - Q ) 2 2/22/3 - a 2 2/32/2)xixi + h-c + ••• , 
where the ellipses denote parts not relevant for our computations. 



(E.l) 



15 To keep the expression as compact as possible we here denote c9+ with dot and d— with prime. 
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